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Abstract 

Providing physical-layer security for mobile users in future broadband wireless networks is of both 
theoretical and practical importance. In this paper, we formulate an analytical framework for resource 
allocation in a downlink OFDMA-based broadband network with coexistence of secure users (SU) 
and normal users (NU). The SU's require secure data transmission at the physical layer while the 
NU's are served with conventional best-effort data traffic. The problem is formulated as joint power 
and subcarrier allocation with the objective of maximizing average aggregate information rate of all 
NU's while maintaining an average secrecy rate for each individual SU under a total transmit power 
constraint for the base station. We solve this problem in an asymptotically optimal manner using dual 
decomposition. Our analysis shows that an SU becomes a candidate competing for a subcarrier only 
if its channel gain on this subcarrier is the largest among all and exceeds the second largest by a 
certain threshold. Furthermore, while the power allocation for NU's follows the conventional water- 
fiUing principle, the power allocation for SU's depends on both its own channel gain and the largest 
channel gain among others. We also develop a suboptimal algorithm to reduce the computational cost. 
Numerical studies are conducted to evaluate the performance of the proposed algorithms in terms of the 
achievable pair of information rate for NU and secrecy rate for SU at different power consumptions. 
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Orthogonal frequency division multiple access (OFDMA), physical-layer security, secrecy rate, dual 
decomposition. 
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I. Introduction 

Security is a crucial issue in wireless systems due to the broadcasting nature of wireless 
radio waves. It also attracts increasing attention because of the growing demand of private 
data transmission such as online transaction and personal medical information. Traditionally, 
cryptography undertakes most of security work on upper layers, which is based on computational 
complexity. In the standard five-layered protocol stack, security approaches are designed on every 
layer except physical layer. Thus, establishing physical-layer security is of both theoretical and 
practical significance. In this study, we aim to provide physical-layer security for mobile users in 
future broadband wireless networks and formulate an analytical strategy for resource allocation 
to achieve this goal. 

Information-theoretic security provides possibility of secure transmission in the physical layer. 
By exploring secrecy capacity and coding technique, messages can be sent without being decoded 
by any eavesdropper. Information-theoretic security originates from Shannon's notion of perfect 
secrecy [HI. He presented the general mathematical structure and properties of secrecy systems. 
The concept of information-theoretic security and wire-tap channel was defined by Wyner ^ 
and later by Csiszar and Komer [3J, who proved the existence of channel coding that makes the 
wiretapper to obtain no information about the transmitted data. Then the study on information- 
theoretic security is extended to various kinds of channels. Leung et al flU focused on Gaussian 
wire-tap channel and showed that secrecy capacity is the difference between the capacities of the 
main and wire-tap channels. Barros et al f5| studied secrecy capacity in slow fading channels and 
introduced outage into secrecy issues for the first time. In |l6l, Li et al investigated independent 
parallel channels and proved that the secrecy capacity of the system is the summation of the 
secrecy rate achieved on each independent channel. More recently, Zhu et al [7\ studied the 
cooperative power control for secret communications by using artificial noise in symmetric 
Gaussian interference channel. 

Orthogonal frequency division multiple access (OFDMA) has evolved as a leading technology 
in future broadband wireless networks, such as 3GPP Long-Term Evolution (LTE) and IEEE 
802.16 WiMAX. It enables efficient transmission of a wide variety of data traffic by optimizing 
of power, subcarrier or bit allocation among different users. In the past decade, tremendous 
research results have been reported on the resource allocation of OFDMA downlink networks. 
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such as [[8l- lfT3l . where the problem formulation differs mostly in optimization objectives and 
constraints. The work in |[8l appeared as one of the earliest results on margin adaptation for 
minimizing the total transmit power with individual user rate requirements. Jang and Lee in 
studied rate adaptation for system sum-rate maximization subject to a total transmit power 
constraint. In ifTOll . Tao et al considered a heterogenous network and studied the power and 
subcarrier allocation problem for maximizing the sum-rate of non-delay-constrained users while 
satisfying the basic rate requirement of each delay-constrained user under a total transmit power 
constraint. Cross-layer optimizations considering utility-function and traffic arrival distribution 
for OFDMA networks were also studied in several works, e.g. [[TT | - [[T3l . Nevertheless, none of 
these works take into account the security issue, which attracts increasing attention recently in 
wireless networks as aforementioned. 

Secrecy or private message exchanges between mobile users and base station (BS) are generally 
needed in present and future wireless systems. Hence, it is essential to consider the security 
demand when assigning radio resources to all users. In lfT4l . the authors made the initial attempt to 
find the power and subcarrier allocation in an OFDM-based broadcast channel with the objective 
of maximizing sum secrecy rate. However, this work is confined to two users only and does not 
consider the coexistence of other types of users. 

In this study, we introduce two types of users according to their secrecy demands in down- 
link OFDMA based broadband networks. The first type of users have physical-layer security 
requirements and should be served at a non-zero secrecy rate. These users are referred to as 
secure users (SU). The other type of users have no confidential messages and do not care about 
security issues. Their traffic is treated in a best-effort way. These users are regarded as normal 
users (NU). All SU's and NU's are legitimate users in the network and have their own data 
transmission with the BS. They are completely honest and always feedback the correct channel 
condition to the BS. Moreover, each of them is equipped with single antenna and only passively 
listens, rater than actively attacking using, for example, multiple antennas or colluding, when 
being a potential eavesdropper of the confidential messages for SU's. The aim of this study is to 
investigate the power and subcarrier allocation problem in such an OFDMA broadband network 
where the BS needs to simultaneously serve multiple SU's and multiple NU's. 

The resource allocation problem in this paper possesses major differences compared with those 
without secrecy constraint, owing to the coexistence of SU's and NU's. Firstly, the legitimate 
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subcarriers to be assigned to an SU can only come from the subcarrier set on which this SU 
has the best channel condition among all the users. This is because for each SU, any other user 
in the same network is a potential eavesdropper. The authors in IfTSl showed that the secrecy 
capacity of a fading channel in the presence of multiple eavesdroppers is the difference between 
the capacities of the main channel and the eavesdropper channel with the largest channel gain 
among all eavesdroppers. As a result, a non-zero secrecy rate on a subcarrier is possible to 
achieve only if the channel gain of the SU on this subcarrier is the largest among all the users. 
Note that this observation is very different from that in conventional OFDMA networks where 
a user is still able to occupy some subcarriers and transmit signals over them even if its channel 
gains on these subcarriers are not the largest. Secondly, even if an SU has the best channel 
condition on a given subcarrier, assigning this subcarrier to the SU may not be the optimal 
solution from the system perspective when the quality of service of NU's is taken into account. 
This is due to the fact that the achievable secrecy rate of channels with Gaussian noise is the 
subtraction of two logarithmic functions [4J, i.e. 

Cs = [log(l + P«m) - log(l + PaE)]"- (1) 

where P is the transmit power, a^i and aE are the channel-to-noise ratios (CNR) of the main 
channel and eavesdropper channel, respectively, and [x]"*" = max{0,x}. In the case of large P, 
if the gap between aM and a^; is not large enough, the achievable secrecy rate can be rather 
low. In this case it may be more beneficial for the system to assign the subcarrier to an NU for 
transmitting non-confidential data traffic. 

Our goal is to find an optimal power and subcarrier allocation policy to maximize the long- 
term aggregate information rate of all NU's while maintaining a target average secrecy rate 
of each individual SU under a total power constraint. Since a non-zero instantaneous secrecy 
rate for each SU cannot be guaranteed all the time due to channel fading, we assume that the 
transmission of confidential messages can wait until the channel condition of the SU becomes 
favorable. As a result, the average secrecy rate requirement instead of instantaneous secrecy rate 
requirement is considered for each SU. In addition, to make our problem more complete, both 
long-term average and peak total transmit power constraints at the BS are considered. 

Finding the optimal power and subcarrier allocation policy with respect to channel conditions 
of all users in the considered system is a functional optimization problem. We solve this 
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problem in dual domain using the decomposition method in an asymptotically optimal manner. 
In particular, the joint subcarrier assignment and power allocation at given dual variables can 
be determined on a per- subcarrier basis. Our analysis shows that for each subcarrier, if the 
user with the largest channel gain belongs to the category of NU's, then this subcarrier will be 
assigned to this NU. Otherwise if it is an SU, then this subcarrier will only be assigned to this 
SU when its channel gain exceeds the second-largest channel gain by a certain threshold. This 
analytical finding agrees with the intuition mentioned above. It is also shown that while the power 
allocation across the subcarriers assigned for NU's follows the standard water-filling principle, 
the power allocation for SU's depends on both the SU's channel gain and the second-largest 
channel gain. Based on the insight derived from the optimal policy, we further propose a low- 
complexity suboptimal power and subcarrier allocation algorithm. Firstly, we assign resources 
to only SU's as if all the NU's were pure eavesdroppers without data traffic. Then, the residual 
resources are allocated to NU's. The underlying mechanism of this algorithm is to decouple the 
joint update of the Lagrangian multipliers as required in the optimal algorithm so as to acquire 
a linear complexity in the numbers of users and subcarriers. 

The rest of the paper is organized as follows. Section II describes the system model and 
problem formulation. In Section III, we introduce dual decomposition method to solve this 
problem under average power constraint. The same problem under peak power constraint is solved 
in Section IV. An efficient suboptimal resource allocation algorithm is proposed in Section V. 
Section VI presents the performance of the proposed scheme via simulation. Finally, we conclude 
this paper in Section VII. 

II. System Model and Problem Formulation 

We consider the downlink of an OFDMA broadband network with one BS and K mobile 
users. The first Ki users, indexed as k = 1,. . .,Ki, are SU's. Each of them has confidential 
messages to communicate with the BS and therefore demands a secrecy rate no lower than a 
constant Ck, for 1 < k < Ki. The other K — Ki users, indexed as k = Ki + 1, . . . , K, are 
NU's and demand service of best-effort traffic. NU's and SU's are all assumed to eavesdrop the 
legitimate channel non-cooperatively and each has only one antenna. The communication link 
between the BS and each user is modeled as a slowly time-varying frequency- selective fading 
channel. The channel coefficients remain approximately unchanged during each time frame, but 
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vary from one frame to another in a random manner. The total bandwidth is logically divided into 
orthogonal subcarriers by using OFDMA with each experiencing slow fading. As a central 
controller, the BS knows the channel information of all users, finds the allocation policy and 
then assigns power and subcarriers to mobile users at each transmission frame according to the 
instantaneous channel state information (CSI) of all users. The total transmit power of the BS 
is subject to either a long-term average or peak constraint, both denoted as power constraint P. 
We assume full statistical knowledge and instantaneous knowledge of CSI at the BS and that 
each subcarrier is occupied only by one user at each time frame to avoid multi-user interference. 

Let „ denote the CNR of user k on subcarrier n for all k and n. The system channel 
condition is denoted by the set a = {afc,„}, which have a joint probability density function 
of f{a). Let Q{a.) = {^li, ...,^Ik} denote the subcarrier assignment policy with respect to 
the system channel condition a, where represents the set of subcarriers assigned to user 
k. Furthermore, let p(a) = V/c, Vn} denote the corresponding power allocation policy, 

where pk,n represents the transmit power allocated to user k on subcarrier n. We next present 
the achievable secrecy rate of SU and the achievable information rate of NU separately, under 
a given resource allocation policy {p(a), i7(a)}. 

For SU k,l < k < Ki, since the subcarriers are parallel to each other, the achievable secrecy 
rate at a given channel realization is the summation of those achieved on each subcarrier in the 
presence of A' — 1 potential eavesdroppers, and can thus be expressed as ifTSl : 

rl=Yl 'In (2) 

where 

'fk,n = [log(l + Pk,nak,n) " log(l + Pfe,n/3fe,n)]^ • (3) 

Here /^^ „ = max ak',n denotes the largest CNR among all the users except user k on subcarrier 
n. The expression ([3]) means that non-zero instantaneous secrecy rate for SU A; on a subcarrier 
is possible to achieve only if its CNR on this subcarrier is the largest among all the K users. 
On the other hand, the achievable information rate of NU k for Ki < k < K is given by: 

rk= ^ rk,n, (4) 

where 

rk,n = log(l + Pk,n(yk,n)- (5) 
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The problem is to find the optimal power and subcarrier allocation policies {p{cx), {Q{cx)} so 
as to maximize the average aggregate information rate of the K — Ki NU's while satisfying the 
individual average secrecy rate requirement for each of the Ki SU's. We consider both peak and 
average power constraints. This functional optimization problem can be expressed as 

max E Wfc rk,n (6) 

subject to E ( ^ „ ) > C^, 1 < A; < i^i (7) 

K 

Yl p^'- ^ p 

k=l nGfifc 

or ^{yY. P^A ^ P (9) 
Pfc,n>0>VA;,n (10) 
fiiU...UfiA' C {1,2, ...,iV} 

VLk are disjoint (11) 

where notation E represents statistical average over the joint distribution of channel conditions, 
i.e. E[-] = / (■)f(a)dcx, and cok is a weighting parameter of NU k, representing its quality-of- 
service demand. Constraint ([8]) is the peak total power constraint while ^ is the average total 
power constraint. 

In the above formulation, we impose the long-term average secrecy rate constraints because 
instantaneous non-zero secrecy rate cannot be guaranteed at every frame. It is possible due to 
channel fading that in a certain frame, an SU is the best user on none of the N subcarriers and 
thus obtains a zero secrecy rate. Therefore, power and subcarrier allocation should be adapted 
every frame to meet a long-term secrecy rate requirement. 

III. Optimal Resource Allocation Under Average Power Constraint 

In this section, we solve the problem with the average power constraint formulated above. It is 
not difficult to observe that this problem satisfies the time-sharing condition introduced in [|T6l . 
That is, the objective function is concave and constraint ([7]) is convex given that rl is concave 
in „ and that the integral preserves concavity. Therefore, similar to the OFDMA networks 
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without secrecy constraint, we can use dual approach for resource allocation and the solution is 
asymptotically optimal for large enough number of subcarriers. 

Define 'P(a) as a set of all possible non-negative power parameters {pk,n} at any given system 
channel condition ct satisfying that for each subcarrier n only one „ is positive. This definition 
takes into account both the power constraint (flOl) and the exclusive subcarrier allocation constraint 
(fTTI) . The Lagrange dual function is thus given by 

" K N 
_k=Ki+l n=l 



max 



^/ife (e 

k=l 



N 



n=l 



X \ P-E 



K N 



,k=l n=l 



[OL] 



(12) 



where = (/ii, . . . , hki) h and A > are the Lagrange multipliers for the constraints (|7]) and 
dll) respectively, and notation E stands for the statistical average over all channel conditions a. 
Then the dual problem of the original problem Q is given by 

min g{fj,,\) 

s.t. /X ^ 0, A > 0. (13) 

Observing (fT2l) . we find that the maximization in the Lagrange dual function can be decom- 
posed into independent sub-functions as: 



N 



n=l 



k=l 



where 



with 



QnifJ., A) = max E [J„ (/x, A, a, {pk,n}k)] 



(14) 



(15) 



K Kl K 

J„(/X, A, a, {pfc,„}fc) = ^ ^fc^fc.n + Atfc^fc.n ~ Pk,n- (16) 

k=Ki+l k=l k=l 

For fixed /j, and A, the maximization problem in (fTSi) is a single-carrier multiple-user power 
allocation problem. Given that the expectation E is over the channel condition a, we can obtain 
the maximum by directly maximizing the function (fT6l) . Based on the above sub-problems, we 
now discuss the optimality conditions of power allocation and subcarrier assignment, respectively 
in the following two subsections. 
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A. Optimality Condition of Power Allocation 

The function in (fT6l ) is concave in „ and hence its maximum value can be found by using 
Karush-Kuhn-Tucker (KKT) conditions. Specifically, suppose that subcarrier n is assigned to user 
k, then taking the partial derivation of A, o:, {pk,n}k) with respect to pk,n and equating it 

to zero, we obtain the following optimality condition of power allocation: 

1 1 \ / 1 



Pk,n 



1 1 

^(^k,n Pk.n 



4/ifc 

A VA 



^k,n 



1 

«fc,n 



1 



C^k,n Pk,n 



for k 



Ki, and 



Pk,n 



UJk_ 

A 



1 



Qik,n 



(17) 



(18) 



for k = Ki + 1,...,K. 

We can conclude from (fTSl) that the optimal power allocation for NU's follows the conventional 
water-filling principle and the water level is determined by both the weight of the NU and 
the average power constraint. On the other hand, it is seen from (flTl) that the optimal power 
allocation for SU's has the same form as the result obtained in [fTTl for conventional fading 
wire-tap channels, as expected. By observing (flTl) closely it is also seen that the SU must satisfy 
«fc n — /3fc n > — in order to be allocated non-zero power. This means that the power allocation 
for SU depends on both the channel gain of the SU and the largest channel gain among all the 



other users. Moreover, it is non-zero only if the former exceeds the latter by the threshold 



B. Optimality Condition of Subcarrier Assignment 

Next, substituting (flTI) and (fTSi) into (fT5l) and comparing all the K possible user assignments 
for each subcarrier n, we obtain 



5(„(/x, A) = E 



where the function Hkn{-) is defined as 



max Hu „ (u, A, a.) 

l<k<K ' ^ 



1 + Plnf^Kn 



(19) 



(20) 



for 1 < A; < and pi „ defined in (fTTl ). and 



iJfc.n (m, a, a) = Uk 



log 



^kak,n 

A 



^k 



(21) 
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for Ki<k<K. 

From (fT9l) it is observed that the function Hk^n defined in (l20l) and (l2T]) plays an important 
role in determining the optimal subcarrier assignment. In specific, for any given dual variables 
fi and A, the subcarrier n will be assigned to the user with the maximum value of iffc,„- That 
is, the optimality condition for subcarrier assignment is given by 



/c* = aigma.xHk n, for n = 1, N. 



(22) 



Note that for k = Ki + 1, K, Hk,n is monotonically increasing in a^.n- Therefore, the NU 
with larger „ is more likely to be assigned subcarrier n. We also notice that for k = 1, Ki, 
Hk,n > only when SU k has the largest a^.n among all the K users and satisfies ak,n > 
(3k,n + Otherwise, Hk^n = 0. In other words, an SU becomes a candidate for subcarrier n 

only if its CNR is the largest and is \/fik larger than the second largest. 



C. Dual Update 

Substituting QnifJ^, A) for n = 1, . . . , N into (fT4l) . we obtain 5f(/x, A). As studied in ifTSl . the 
dual problem (fT3l) is always convex and can be minimized by simultaneously updating (/x, A) 
using gradient descent algorithms. However, note that g(fj,,\) is not differentiable due to the 
discontinuity of subcarrier assignment and hence its gradient does not exist. Nevertheless, we 
present a subgradient of g{n, A) as follows: 



E 



N 



n=l 



K N 



A\ = P-E 

_k=l n=l 

where r|*„ is obtained by substituting the optimal pi „ into ([3]). A brief proof is as follows 
From the expression of g{fi, \) in (fT2l) . we get 



(23) 
(24) 



^?(m',a') 



> E 



K 



N 



,k=Ki+l n=l 



k'n 



k=l 



N 



n=l 



Ck 



X' \ P-E 



K N 



fc=i 



N 



n.=l 



(A' - A) P - E 



,fc=l n=l 
,k=l n=l 



.(25) 



The results in (1231) and (124)) are thus proved by the definition of subgradient. 
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In the case when the numerical computation of statistical average is too complex, we can 
change to time average when the channel fading process is ergodic. Specifically, 



After finding the optimal dual variables {/i^} and A*, the optimal power and subcarrier 
allocation policy is then obtained by substituting them into the optimality conditions (fTTI) . (fTSi) 
and (l22l) . 

D. Discussion of Feasibility 

To make the above optimization problem feasible, the secrecy requirement Ck for each SU 
must be chosen properly according to the total power constraint P. In this subsection, we will 
derive an upper bound of average secrecy rate each SU can obtain. If Ck is set equal or greater 
than this upper bound, the secrecy requirement cannot be satisfied and thus this optimization 
problem is not feasible. 

For simplicity, we assume that the channel conditions of all the K users on each subcarrier 
are independently and identically distributed and follow Rayleigh distribution. Then, each SU 
has a probability of \/K to be the best user on each subcarrier. As a result, we can write the 
average achievable secrecy rate of each SU k as 



Assuming that the transmit power goes to infinity, the maximum per- subcarrier secrecy rate 
can be obtained as 




(26) 




(27) 



Rl = -E[rlJ. 



(28) 



1 + 



lim r 



k,n 



log 



(29) 
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Therefore, an upper bound on Rl can be theoretically derived using order statistics as 

1^1 



OO POO 

/ log—f{vi,V2)dvidv2 

Jv2 ^2 

N r°° r°° _-2 N~2i ^2 1 -1 

— / / A^(iV-l)(l-e p) -e p -e p \og—dvidv2 

K Jo Ju2 P P ^2 

poo poo 

K 



N 
K 





OO POO 



Af-2 1 1 , 

e p j -e p -e p \0gu1dv1du2 
P P 



N-2\ 1 

A^(A^— 1)(1 — e p) -e p -e p \ogiJ2di'idv2- 



(30) 



'0 JV2 

In the above derivation, ui and 1^2 denote the largest and second largest CNR's on a given 
subcarrier, /(z/i, h'2) is the joint probability density function of ui and U2^ which can be obtained 
through order statistics. We assume Rayleigh fading and the probability distribution function of 



CNR is /(x) 



-e p 



, where p is the mean value. Since it is difficult to further express the 



integrals in (1301) in a closed form, we use numerical integration to get the upper limit value of 
secrecy rate. As an example, when the parameters are set to be = 64, K = 8, Ki = 4 and 
p = 1, we obtain Rl < 3.5nat/OFDMA symbol. In this case, when the secrecy rate constraint 
Cfc > 3.5, the problem becomes infeasible. 



IV. Optimal Resource Allocation under Peak Power Constraint 

In Section III, we solved the optimization problem under long-term average power constraint. 
In this section, we solve the same problem Q but under peak power constraint. Note that the 
peak power constraint is more suitable for practical systems as the BS usually has a maximum 
radiation power. By using the similar definition of power parameter set V{cx) as in the previous 
section, the associated Lagrange dual function is given by 



max 

{Pk,n}eV{a) 



K 



N 



+ 



=A'i + l n=l 
V k=l n=l J J 




N 



n=l 



c. 



(31) 



The difference from the dual function with average power constraint in (1121) is that the dual 
variable A(q;) associated with the power constraint ([91) is a function of the system channel 
condition a. 
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This dual function can be similarly decomposed into independent sub-functions with each 
given by 

gnifJ-, X{cx)) = niax E[J„(/x, A(q:),q:, {pfc,n}fe)]> (32) 

{Pfc,n}6P(a) 

with 

K Ki K 

Jnif^,Xia),CX,{pk,n}k) = Wfcrfc,n + ^/ifcr^_„ - A(q:) ^Pfc,n- (33) 

k=Ki+l k=l k=l 

Given that the order of expectation operation and the max operation can be reversed, we can 
obtain 5'„(/^, A (a)) by maximizing the function (l33l) for every ak,n- 

Suppose that subcarrier n is assigned to user k. Taking the partial derivative of J„(a*5 A(q;), ol, {pk,n}k) 
with respect to „ and equating it to zero, we obtain the optimality conditions of power 
allocation 



for k 




Pk,n 



4/ifc 



k,n 



Pk,n 



UJk 



(34) 



(35) 



_X{a.) ak,n_ 
for k = Ki + l,...,K. 

Comparing (|34l) and (|35l) with (flTI ) and (fTSi) . we observe that the optimal power allocations 
under the average power constraint and the peak power constraint have similar structure. The 
difference lies in the dual variable that controls total power. For average power constraint, A is 
a constant for all ak,n- For peak power constraint, this variable changes with {afc,n}- 

The rest of the algorithm is alike to the problem in the previous section. The difference is 
that the subgradient to update A (a) is changed to 



K N 



k=l n=l 



(36) 
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V. SuBOPTiMAL Power and Subcarrier Allocation Algorithm 

The complexity of the optimal power and subcarrier allocation policy presented in the previous 
two sections mainly lies in the joint optimization of Lagrange multipliers /x = (/ii, . . . , ^Ki) 
and A. If we choose the ellipsoid method, it converges in 0{{Ki + l)^logi) iterations where 
e is the accuracy ifTSl . Thus, if the number of SU is large, the computational complexity may 
not be favorable for practical implementation. Based on the insight derived from the optimal 
power and subcarrier allocation policy in the previous sections, we present in this section a low- 
complexity and efficient suboptimal power and subcarrier allocation algorithm. For simplicity, 
only the average power constraint is considered. The extension to the peak power constraint 
is simple. The idea of this scheme is to first assign the resources to only SU's as if all the 
NU's were pure eavesdroppers without data transmission. After that, the residual subcarriers 
and power, if any, are distributed among NU's. By doing this, the joint update of the Lagrange 
multipliers will be decoupled as detailed below. 

In this scheme, the power allocation adopts the expressions in (flTI) and (fTSi) except that the 
parameter \/nk,k = 1, . . . , Ki in (flTl) is replaced by a new variable z/^. Also, in (fTSi) we define 
Lfc = WfcLo, for k = Ki + 1, . . . , K) where Lq = \. The power allocation scheme among NU's 
is water-filling with different water levels which are proportional to NU's weights. Through 
simple observation, for SU pk,n is monotonically decreasing in Uk and thus r|, „ is monotonically 
decreasing in Uk- Intuitively, for NU „ is monotonically increasing in water level Lq. Therefore, 
Uk and Lq can be found through two separate binary searches. We use a set of training channel 
realizations to compute the statistical average numerically. The set is sufficiently large to assure 
that their distribution converges to the statistical distribution. 

The outline of this suboptimal algorithm is presented below. 

Suboptimal Algorithm 
Find the optimal Vk to achieve the secrecy rate requirement for /c = 1, . . . , Ki. 

1) Set v}!^ sufficiently large, vj^^ = and Uk = \{vk^ + i']^^). 

2) For every training channel realization a. 

Find Vtk = {n : ak,n > maxfc/^fc(afc/,„) + Uk}; 

Compute pk{a.) = Yl Pfc,n(«) and rl{a) = Yl i"l{a), where Pfe,„(a) and r^._„(a) are 
computed according to (fTTI) and (HI), respectively with X/nk replaced by z/^; 
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Compute r| = E [r^(a)] and = E [pk{a.)]. 

3) If rl > a, u^"" = Uk, else uj^^ = u^. Set z/,. = + uj^^) 

4) Repeat Steps 2)-3) until \fl - Ck\ < eCk, for each k e [1, Ki] 

— _ 

5) Compute the power consumed by SU's, Psu = Yl Pk- 

k=l 



Find the optimal water level Lq to meet the power constraint. 



6) Set L^^ sufficiently large, L^^ = and Lq = + L^^ 



7) For every training channel realization ol 

For every residual subcarrier n ^ [j Q^, i.e. not occupied by SU's 

k=l 

Find k = arg max H^n according to (|2T]) with 1/A replaced by Ln; 

k£{Ki,K] 

For the found k, compute Pk,n{oi) and rfc ,^(a) according to (fTSi) and Q with 1/A 



replaced by Lq; 



Compute Tfc = E 



AT 



E ^/fc,n(Q:) 



and Pnu = IE 



ir N 

E E Pfc,n(Q:) 

k=Ki+l n=l 

8) If Pr,u>P- Psu, Lll^ = Lo, else L^^ = Lq. Set Lq = + 



n=l 



9) Repeat Steps 7)-8) until \Psu + Pnu - P\ < ^P- 



In the algorithm, we first find the the subcarrier set ^k{k = 1, . . . , Ki) of SU's. The criterion 
is whether SU's CNR is larger than the second largest. As fi^'s are disjoint, the optimal Uk 
satisfying Ck can be obtained through Ki independent binary searches. After getting {I'k}, the 
power allocation for the Ki SU's is determined and the total power left for K — K\ NU's is 
also known. The water levels of power allocation (fTSi) for NU's can be searched until the rest 
of the total power for NU's is used up. 

Since v^ik = 1, . . . , Ki) and Lq are obtained individually by binary search, this suboptimal 
algorithm converges in 0{Ki log ^ + log ^) = 0{{Ki + 1) log ^) iterations. Note that the optimal 
algorithm converges in 0{{Ki + l)^logi) iterations as mentioned in the beginning of this 
section. In addition, the computational loads in each iteration of both the optimal and suboptimal 
schemes are linear in KN\a.\ where |a| is the number of the training channel realizations. So 
the suboptimal scheme reduces the complexity by about 7^^- 



VI. Numerical Results 

In this section, we provide some numerical results to evaluate the performance of proposed 
optimal and suboptimal resource allocation algorithms under both average and peak power 
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constraints. In the simulation setup, we consider an OFDMA network with = 64 subcarriers 
and K = 8 mobile users, among which Ki = 4 are SU's and K—Ki = 4 are NU's. For simplicity, 
all the weighting parameters cj^'s for NU's are set to 1 and the secrecy rate requirements C^'s 
for SU's are set to be identical, denoted as Ck = Rsu- Let R^u denote the average total 
information rate of the NU's. Here we use nat instead of bit as the measurement unit of data 
for computational convenience. The channel on each subcarrier for each user is assumed to 
be independent and identically distributed Rayleigh fading with unit mean-square value for 
illustration purpose only. Note that more sophisticated multi-path broadband channel models, 
such as HiperLan/2 channel model A, can be easily applied since our analytical framework is 
general and applicable to arbitrary channel distributions. The system total transmit SNR defined 
in the simulation is equivalent to the average power budget P (or peak power budget if the peak 
power constraint is concerned) on all the subcarriers at the base station, assuming unit noise 
power. 

To evaluate the optimal and suboptimal power and subcarrier adaptation schemes, we introduce 
two non-adaptive schemes in this simulation as benchmarks. In these two non-adaptive schemes, 
subcarrier assignment is fixed beforehand while power allocated to the predetermined subcarrier 
sets conforms to (flTl) and (fTSi) . In the first fixed subcarrier assignment scheme, denoted as FSA- 
1, the 64 subcarriers are equally assigned to the 8 users regardless of user type and thus each 
user obtains 8 subcarriers. In the second scheme, denoted as FSA-2, the SU's are given higher 
priority and each is assigned 12 subcarriers, whereas each NU is assigned 4 subcarriers. 

We first demonstrate the pair of achievable average total information rate of NU's and feasible 
average secrecy rate requirement of each individual SU {Rnu,Rsu) at a given average total 
power constraint. Fig. [T] shows the rate pair at fixed total transmit SNR=30dB. First it is observed 
that using both optimal and suboptimal algorithms, R^u decreases with the increase of secrecy 
requirement Rsu and falls sharply to zero at around Rsu = 3.5 nat/OFDM symbol. Recall the 
discussion in Section III-D where an upper bound on the secrecy rate requirement to make the 
problem feasible was obtained as 3.5 nat/OFDM symbol in the case of = 64 subcarriers and 
K = 8 users. This explains the drastic fall of Rnu- It is then observed from Fig. \T\ that the 
suboptimal algorithm only incurs less than 20% loss in Rj^u when achieving the same Rsu 
compared with the optimal algorithm. Now comparing the optimal and suboptimal schemes 
with the non-adaptive ones, both of them earn great advantage over the two benchmarks FSA-1 
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and FSA-2. In particular, the maximum feasible points of the two fixed subcarrier assignment 
methods appear at around Rsu = 0.44 and 0.66 nat/OFDM symbol, respectively, which are 
much lower than those in the optimal and suboptimal algorithms. 

Fig. [2] and Fig. [3] show, respectively, the average power consumption and the average number 
of subcarriers assigned to all SU's with respect to different Rsu when the total average transmit 
power is fixed as 30 dB. From Fig. [2l we notice that the optimal scheme spends more power on 
SU's than the suboptimal one. It is observed from Fig. [3] that the number of occupied subcarriers 
by SU's increases with the growing of Rsu and reaches 32 at the feasible point of Rsu = 3.5 
nat/OFDM symbol for both the optimal and suboptimal schemes. Note that 32 is an expected 
number because when Rsu is very close to the feasible point, adaptive schemes tend to assign as 
many subcarriers as possible to SU's and on average, all SU's can occupy half of total subcarriers 
at most. Additionally, the optimal scheme assigns less subcarriers to SU's than the suboptimal 
one. For FSA-1 and FSA-2, the number of occupied subcarriers by SU's is also increasing with 
the feasible Rsu and is smaller than the number of pre-assigned subcarriers to SU's, i.e. 8 and 
12, respectively for FSA-1 and FSA-2. This indicates that fixed subcarrier assignments waste 
subcarriers compared with adaptive ones. Note that the number of subcarriers assigned to all 
NU's for the optimal and suboptimal schemes can be straightforwardly obtained by subtracting 
the numbers shown in Fig. [3] from the total number of subcarriers, 64. 

We next demonstrate the relation between R^u and total transmit SNR for a given Rsu =0.4 
nat/OFDM symbol in Fig. H Note that the constraint Rsu = 0.4 is feasible when SNR> -2dB 
for the optimal and suboptimal schemes and when SNR> 3 or 9 dB, respectively for FSA-1 
and FSA-2. From Fig. IH we first observe that at low SNR region (SNR<18dB), the optimal 
and suboptimal schemes perform close to each other. As SNR becomes larger, the suboptimal 
scheme only incurs a marginal performance loss. It is also seen that the function curves of FSA- 
1 and FSA-2 intersect at about 14dB. When the transmit SNR is lower than 14dB, FSA-2 is 
superior to FSA-1 in terms of R^u- This is because when power is limited, FSA-2 assigns SU's 
more subcarriers and thus save more power for NU's. If transmit SNR is higher than 14dB, 
FSA-2 becomes inferior to FSA-1. The reason lies in that when power is sufficient to meet 
secrecy rate requirements, FSA-2 wastes subcarriers on SU's and leaves NU's less to promote 
their information rates. 

Finally, in Fig. [51 we compare the performance under average and peak power constraints and 
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show the achieved aggregate rate of NU's at different SNR. It is observed that under the two 
power constraints, the two curves differ slightly in low SNR region and almost coincide at high 
SNR region for both optimal and suboptimal schemes. 

To conclude the above results, the proposed optimal and suboptimal resource allocation 
schemes significantly outperform those with fixed subcarrier assignment. This observation in- 
dicates the great importance of carefully coordinating the subcarrier allocation with adaptation 
to the system channel conditions. Moreover, the suboptimal scheme provides a good tradeoff 
between performance and complexity. 

VIL Conclusions and Discussions 

This work investigate the power and subcarrier allocation policy for OFDMA broadband 
networks where both SU's and NU's coexist. We formulate the problem as maximizing the 
average aggregate information rate of NU's while satisfying the basic average secrecy rate 
requirements of SU's under either an average or a peak transmit power constraint. We solved the 
problem asymptotically in dual domain by using decomposition method. Results show that the 
optimal power allocation for an SU depends on both its channel gain and the largest channel gain 
among others. We also observe that an SU becomes a valid candidate competing for a subcarrier 
only if its CNR on this subcarrier is the largest among all and larger enough than the second 
largest CNR. To reduce the computational cost, a suboptimal scheme with favorable performance 
is presented. Numerical results show that the optimal power and subcarrier allocation algorithm 
effectively boosts the average total information rate of NU's while meeting the basic secrecy 
rate requirements of SU's. It is also shown that whether it is peak or average power constraint, 
the system performance does not differ much given sufficient number of subcarriers used. 

Before finishing the paper, we provide some discussions. Firstly, we assumed throughout this 
work that the channel state information of each user obtained at the BS is accurate. If a user 
deliberately lies and reports a lower channel-to-noise ratio on certain subcarriers, it would get 
higher chance of eavesdropping SU's private message and therefore cause secrecy rate loss. 
However, on the other side, its own average information rate or secrecy rate would also be 
reduced due to the less assigned radio resources. Hence, we argue that there is no incentive for 
the users to lie about their channel condition. If, however, the network contains purely malicious 
eavesdropper that does not care about its own transmission, the security in the network can be 
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circumvented by sending the BS false channel measurements. Secondly, if eavesdroppers can 
collude and exchange outputs to decode the message, the network can be regarded as a single- 
input multiple-output system, where there is only a single eavesdropper with multiple receive 
antennas. In this case, our algorithms can still apply except the change that the secrecy rate of an 
SU depends on the sum of the channel gains of all the eavesdroppers rather than the largest one. 
Furthermore, in the case where a user is equipped with multiple receive antennas, an interesting 
topic for future investigation is to allow the user to report only one antenna to the BS and use 
the remaining antennas to eavesdrop. Lastly, since the proposed resource allocation algorithm 
is a centralized one, full knowledge of CSI is required as commonly assumed in the literature. 
Taking into account practicality and system overhead, we will investigate distributed algorithms 
with partial or local channel knowledge in our future work. 
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Fig. 1. Achievable {Rsu, Rnu) pair at total transmit SNR of 30dB. 



21 



700 



600 



b 500 



-° 400 

"D 

E 

« 

c 300 



0) 

° 200 



100 



— Optimal 
— e — Suboptimal 

FSA-1 

■ FSA-2 



a6'646 ' 



6 




Fig. 2. Average power consumption by all SU's versus i?5i7 at total transmit SNR of 30dB. 
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Fig. 4. Rnu versus total transmit SNR at Rsu = 0.4nats/OFDM symbol. 
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Fig. 5. Rnu versus total transmit SNR at Rsu ~ 0.4nat/OFDM symbol under both average and peak power constraints. 



